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ABSTRACT 
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High-resolution spectra of the hot white dwarf G191-B2B, covering the wave- 
length region 905-1 187 A, were obtained with the Far Ultraviolet Spectroscopic 
Explorer {FUSE). This data was used in conjunction with existing high-resolution 
Hubble Space Telescope STIS observations to evaluate the total H I, D I, O I and 
N I column densities along the line of sight. Previous determinations of A r (D i) 
based upon GHRS and STIS observations were controversial due to the saturated 
strength of the D I Lyman a line. In the present analysis the column density of 
D i has been measured using only the unsaturated Lyman /3 and Lyman 7 lines 
observed by FUSE. A careful inspection of possible systematic uncertainties tied 
to the modeling of the stellar continuum or to the uncertainties in the FUSE in- 
strumental characteristics has been performed. The column densities derived are: 
log N(D I) = 13.40±0.07, loglV(0 1 ) = 14.86i0.07, and log JV(N 1 ) = 13.87±0.07 
quoted with 2 a uncertainties. 

The measurement of the H I column density by profile fitting of the Lyman a 
line has been found to be unsecure. If additional weak hot interstellar components 
are added to the three detected clouds along the line of sight, the H I column 
density can be reduced quite significantly, even though the signal-to-noise ratio 
and spectral resolution at Lyman a are excellent. The new estimate of N { H i) 
toward G191-B2B reads: log./V(H 1 ) = 18.18 ±0.18 (2a uncertainty), so that 
the average (D/H) ratio on the line of sight is: (D/H)= 1.661^6 x 10 -5 (2cr 
uncertainty). 

Subject headings: ISM: abundances : ISM — ultraviolet 


1. Introduction 

Deuterium is thought to be produced in significant amount only during primordial Big 
Bang nucleosynthesis (BBN), and to be thoroughly destroyed in stellar interiors. Deuterium 
is thus a key element in cosmology and in Galactic chemical evolution (see e.g. Audouze 
& Tinsley 1976; Gautier & Owen 1983; Vidal-Madjar & Gry 1984; Boesgaard & Steigman 
1985; Olive et al. 1990; Pagel 1992; Vangioni-Flam h Casse 1994; Vangioni-Flam et al. 
1995; Prantzos 1996; Scully et al. 1997; Casse & Vangioni-Flam 1998; Tosi et al. 1998). 
Indeed, its primordial abundance is the best tracer of the baryonic density parameter of the 


'This work is based on data obtained for the Guaranteed Time Team by the NASA-CNES-CSA FUSE 
mission operated by the Johns Hopkins University. 
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Universe, Qb, and the dec ease in its abundance during galactic evolution traces the cosmic 
star formation rate at various epochs. 

The first (indirect) measurement of the deuterium abundance was carried out using 3 He 
in the solar wind, giving the presolar value D/H~ 2.5±1.0xl0 -5 (Geiss & Reeves 1972). The 
first measurements of the D/H ratio in the interstellar medium (ISM) were reported shortly 
thereafter by Rogerson k York (1973), and their value (D/H)~ 1.4±0.2 x 10 5 has remained 
a landmark average value for the interstellar D/H ratio. Finally direct measurements of the 
primordial (D/H) ratio in low-metallicity material at high redshift have been successfully 
carried out these past few years (e.g., Buries 2001 for a review, and references therein). The 
values derived cluster around (D/H)~ 3 x 10 -5 although with significant dispersion, which 
may or may not be real. Quite similarly the measurements of the (D/H) ratio in the Galactic 
ISM towards hot stars wit h the Copernicus satellite lead to many evaluations of D/H (see 
e.g. York and Rogerson 1976; Vidal— Madjar et al. 1977; Laurent et al. 1979; Ferlet et al. 
1980; York 1983; Allen et al. 1992) which also show dispersion around the above York k 
Rogerson (1973) value. This dispersion has been recently confirmed by IMAPS observations 
(Jenkins et al. 1999; Sonreborn et al. 2000), indicating that the D/H ratio may vary by a 
factor ~ 3 in the solar neighborhood, i.e., within a few hundred parsecs. 

In this paper we present a new determination of the D/H ratio on the line of sight 
to the nearby DA white dwarf (WD) G191-B2B based on observations obtained with the 
Far Ultraviolet Spectroscopic Explorer (FUSE, Moos et al. 2000; Sahnow et al. 2000). This 
paper is one of a series in this volume describing the first results of the FUSE (D/H) program 
in the Local ISM (LISM). This program and its results are summarized in the overview paper 
by Moos et al. (2002). 

Observing white dwarfs has many advantages over hot and cool stars for studying the 
D/H ratio, as explained in Yidal-Madjar et al. (1998): these targets can be chosen close to 
the Sun, in order to avoid a complex line of sight structure, and in the high temperature 
range, so that the interstellar absorption is superimposed on a smooth stellar continuum. 
The risk of contamination by low column density H I fluffs possibly present in the hot star 
winds (Gry, Lamers k Vidal-Madjar 1984) is negligible for WDs, and their hot continuum 
offers the possibility of observing the numerous UV lines of N I and especially O I, which 
is a very useful tracer of 1 1 1 and D I. The (D/H) ratio has already been measured toward 
four white dwarfs, using the HST: G191-B2B (Lemoine et al. 1996; Vidal-Madjar et al. 
1998; Sahu et al. 1999). HZ43A (Landsman et al. 1996), Sirius B (Hebrard et al. 1999) and 
Feige 24 (Vennes et al. 2H00). For HZ43A, Sirius B and Feige 24 the average D/H values 
obtained are compatible v ith the local ISM (LISM) D/H determination (Linsky 1998) made 
within the Local Interstellar Cloud (the LIC in which the sun is embedded), although in the 
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case of Sirius B this compatibility is marginal. 

In the case of G191-B2B, it was found that the line of sight comprises one neutral 
region corresponding to the LIC, and two more ionized absorbing components (Lemoine 
et al. 1996; Vidal-Madjar et al. 1998). The average (D/H) ratio (defined as the ratio of 
the total column densities of D I and H i) was found to be (D/H)= 1.12 ± 0.08 x 10 -5 
(Vidal-Madjar et al. 1998), significantly different from the value measured toward Capella 
(D/H)lic = 1-5±0.1 x 10 -5 (Linsky 1998). The (D/H) ratio measured toward G191-B2B has 
been contested by Sahu et al. (1999), who used STIS data and concluded to the presence of 
two interstellar components only, and to a D/H ratio compatible with that observed toward 
Capella. The disagreement resides in the evaluation of the total D I column density (Vidal- 
Madjar 2000; Sahu 2000), and arises presumably because the Lyman a D line is saturated 
and the column density is thus sensitive to the line profile. The number of components 
assumed on the line of sight may also introduce differences between the analyses of these 
groups (see Vidal-Madjar 2001 for a detailed discussion). 

In the present work we re-examine these issues, making use of high quality FUSE and 
STIS observations of G191-B2B. We first measure the N(D i), N(0 i) and N( N i) column 
densities using the unsaturated lines of these elements in the FUSE datasets, notably Ly- 
man (3 and Lyman 7 for D I (Section 2). We then analyze the recent high quality STIS 
observations in Section 3, and provide explicit evidence for the presence of three absorbing 
components (at least) on the line of sight. We also provide a refined estimate of the total 
H I column density. All throughout this work, considerable effort has been put on quan- 
tifying possible systematic uncertainties related to fixed-pattern noise, detector artifacts, 
background uncertainties, wavelength calibration and modeling of the stellar continua, as 
well as to the velocity structure of the line of sight. In particular, we argue in Section 3.3 
that previous estimates of the total N( H 1) are subject to a large systematic uncertainty 
related to the possible presence of additional weak [iV(H i)< 10 14 cm~ 2 * * ] hot ( T ~ 10 5 K) 
components. This effect may have a large impact on our understanding of the observed 
variations of the (D/H) ratio in the ISM, as it may affect other lines of sight, and is the 
subject of a companion paper (Vidal-Madjar & Ferlet 2002). We provide a summary of our 
results and a short discussion in Section 4; an overall discussion of the FUSE results is given 
by Moos et al. (2002). 


2. Column densities of D I, O I and N 1 measured with FUSE 

The FUSE instrument (Moos et al. 2000; Sahnow et al. 2000) gives access to the 

wavelength range 905-1 187A on eight detector segments, through three different entrance 
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apertures (the large, medium and narrow apertures - LWRS, MDRS and HIRS, respectively) 
and four channels (two LiF and two SiC). In particular, it gives access to the higher order 
unsaturated Lyman lines »f D I and can thus offer a reliable estimate of N(D i). Similarly, 
many weak lines of 0 I and N I are present in this bandpass, as shown in Table 1. The FUSE 
spectrum of G191-B2B in the SiClB channel observed with the HIRS aperture is shown in 
Fig. 1, in w r hich the lines of D I, N I, 0 I used in our study are indicated. 

The white dwarf G1ET-B2B has been observed several times with FUSE through the 
three different entrance apertures; the observation log is given in Table 2. Since the target is 
point-like, the resolutions obtained with these apertures should not be too different, but there 
is a clear variation due to some off-axis effect when moving from one aperture to another. 
The main problem comes from the geocoronal Lyman /? emission whose strength increases 
with the surface area of the aperture, and which blends with the interstellar H I absorption. 
This leads us to employ different approaches in the analysis of the data, depending upon the 
slit used, as discussed below. 

The observations were treated through the version 1.8.7 of the CALFUSE pipeline. The 
data were then collected channel by channel by series of subexposures. Due to the sensitivity 
of FUSE and the brightness of the source, the data were obtained in the HISTOGRAM mode 
(a spectral imaging mode, in which a two-dimensional spectrum was accumulated onboard). 
Each subexposure has a $/N ratio that is high enough to clearly see the strong interstellar 
or photospheric features. The different subexposures could then be easily aligned on top of 
each other to compensate for the slow thermal drifts that displace the wavelength scale from 
one subexposure to another. These drifts never translates into more than ±5 pixels. A high 
S/N spectrum is thus reconstructed in each channel and for each aperture. This process 
preserves the ultimate FUSE resolution and it also partly eliminates some of the detector’s 
fixed pattern noise by act ng like a random FP-SPLIT procedure. 

Since the FUSE sampling is ~ 10 pixels per resolution element, we rebinned the spectra 
by three. We measure column densities by profile fitting of the various lines observed, using 
the code Owens . f developed by one of us (M. L.) and the French FUSE team. Our general 
procedure for data analysis is as follows. We split each spectrum into a series of small sub- 
spectra centered on absorption lines to be analyzed, whose width is a few A depending on the 
line density and fit all lines contained in all sub-spectra simultaneously. Between different 
sub-spectra, the zero-point velocity offset is left to vary to compensate for the wavelength 
calibration ripple effects in FUSE data; the Line Spread Function is generally taken to 
be a Gaussian v r hose FW HM varies freely during the fitting procedure and independently 
from sub-spectrum to sul -spectrum. It is also checked that considering a double Gaussian 
LSF with free wudths, free amplitude ratio and zero pixel offset between the two Gaussian 
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components, does not change the results. The background levels used for each line are those 
evaluated at the bottom of the closest Lyman line (see Fig. 2 for an example), which leads 
to an upper limit for the derived column densities; fixing the background level at the 0 value 
gives access to the lower limit of the evaluated column densities. This uncertainty is taken 
into account in our procedure; however FUSE is operating in first order grating mounts 
and consequently the stray light level is generally quite low, as can be seen in Fig. 2 for 
Lyman 7. This analysis procedure allows to reduce systematic uncertainties tied to fixed 
pattern noise effects or detector artifacts and uncertain calibration, since different lines of 
the same elements observed through different apertures and in different channels are fitted 
simultaneously. For instance, if a given line is subject to an artifact in the dataset, the x 2 of 
the fit in the region around this line, or equivalently the quality of the fit of this particular 
line should stand out in the overall fit if the number of lines analyzed is sufficiently large, as 
is the case here. Moreover, by letting the unknown instrumental parameters vary freely, i.e., 
the LSF shape and width, the background flux, the wavelength zero-point and the shape of 
the continuum, and be optimized simultaneously with the physical parameters that define 
the absorbers, one marginalizes the final result over these unknowns. Therefore the errors 
attached to the uncertainty on these above characteristics are included in the final errors 
quoted for the column densities. Finally, all lines are also analyzed one by one, i.e. each 
independently of all others, in order to check for overall consistency. A detailed discussion 
of this procedure and survey of possible systematic errors tied to FUSE data is given in a 
companion paper by Hebrard et al. (2002). 

2.1. O I and N I 

Since FUSE cannot resolve the velocity structure of the line of sight (see Section 3), 
we have used only unsaturated lines of D I, N I and O I and assumed only one interstellar 
component to be present. Therefore, all measured column densities in this Section are 
integrated over the line of sight. It has been verified that using three interstellar components 
in the fit of the FUSE did not change the conclusions as long as only non-saturated lines 
were considered; an example of single vs multiple component fits is shown in Fig. 2. 

Note that HST observations cannot give a reliable measurement of the O 1 column 
density without a complete knowledge of the velocity and broadening structure of the line 
of sight since the only absorption line available at 1302A is strongly saturated. In contrast 
the N I triplet at 1200A in the GHRS or STIS bandpass is only moderately saturated and 
may provide a precise estimate of the total JV(N i). The various measurements of iV(N 1) 
and N ( O 1) are listed in Table 3. 
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Since the contamination of interstellar lines by weak photospheric lines is a likely pos- 
sibility in the FUSE spectral range due to the high density of atomic lines in this bandpass, 
we excluded from the fit all lines that presented an obviously excessive equivalent width or 
column density relative to the others. As an example as shown on Figure 3 the N I line at 
954.1042A is very probably blended with a photospheric line which shows up at nearly the 
same wavelength in the Nl.TE spectrum calculated with the TLUSTY/SYNSPEC codes (de- 
scribed in Section 3). Such blending becomes particularly important in the case of G191-B2B 
because the total O I or IN I column densities are small, and their absorptions are relatively 
weak features which may be more easily blended with undetected photospheric ones. This is 
much less important for the other FUSE (D/H) targets (see the synthesis of FUSE results by 
Moos et al. 2002, and the individual line of sights studies by Friedman et al. 2002, Hebrard 
et al. 2002, Lehner et al. 2002, Sonneborn et al. 2002 and Wood et al. 2002) which present 
larger column densities, except for the white dwarf HZ43A (Kruk et al. 2002) which presents 
a pure H atmosphere (Barstow, Holberg & Koester 1995; Dupuis et al. 1998) and thus no 
possible blend with photospheric lines. However note that this line selection process may 
induce in the case of G191-B2B a slight underestimation of the total O I or N I column 
densities. 

The O I and N I lines used in this study are listed in Table 1. Some of them were 
rejected either because they are saturated (marked ‘‘strong” in Table 1), or because of a 
blend with a photospheric feature (marked “blend”), or because they are too weak and 
subject to systematic effects of the noise (marked “weak”). Some were kept in one channel 
but not in an other if a nearby detector defect was identified. For this reason, having the 
same spectrum observed in different independent channels was extremely helpful. 

By following the general method depicted above, the final column density determinations 

are : 


log A r (0 l) tot = 14.86 (±0.05) (±0.04) 
log iV(N i) tot = 13.87 (±0.05) (±0.05) 

where the errors are first statistical (2cr) and second systematic (95.5% c.l.). The statistical 
errors were evaluated through the Ay 2 technique, an example of which is given for the iV(D 1 ) 
determination to follow. However, one important comment is of order at this point. The 
total x 2 f° r 311 spectral hues fitted simultaneously (H I, D I, O I and N 1 ) in the case of 
the HIRS observations is equal to 3211.44 for 1843 degrees of freedom (d.o.f.), even though 
the fit appears very satisfactory by eye. This % 2 value lies well above the upper 95.5% 
confidence level limit for statistical fluctuation at 1966.71. We interpret this large x 2 value 
as a result of a systematic underestimate of the magnitude of the individual pixel errors by the 
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instrument. In order to quantify this uncertainty, we have completed low order polynomial 
fits over relatively smooth and fiat continuum sections close to each of the spectral lines 
analyzed. These simple fits gave a x 2 of 688 for 500 degrees of freedom, again above the 
95.5% c.l. upper limit at 563.9. This indicates that indeed the individual pixel errors are 
underestimated by instrument, probably due to the presence of fixed pattern noise in the 
data. The discrepancy in terms of ratio of measured to expected x 2 is not as strong as for 
the best-fit x 2 values, but this may be due to differences in the error levels for pixels located 
on the continuum and those located at the bottom of the saturated Lyman lines. 

In order to compensate for these effects, we have thus decided to rescale all x 2 measured, 
which amounts to an overall rescaling of the individual pixel error estimates. In particular, 
we have chosen to divide all x 2 by 1.84 = 3211.44/1741.64, where 1741.64 corresponds to 
the lowest expected x 2 a t the 95.5% confidence level for 1843 d.o.f. This rescaling rescales 
(increases) the error bars on the physical parameters that we derive since the Ax 2 contour 
levels are also rescaled. In effect, in order to derive the 3a error bar around a single parameter 
minimum value, we search the values of this parameter that lead to an increase in (non- 
rescaled) x 2 of 9 x 1.84 = 16.6 instead of 9. Our choice of rescaling is thus very conservative 
with respect to the final error bars measured. All presented FUSE statistical errors will be 
evaluated following this method. 

The systematic errors quoted above for O I and N I reflect the range of values obtained 
for the best-fit solutions for different apertures, given in Table 3. Note that in this table, the 
errors quoted include the satistical error plus part of the above systematic errors, since it 
accounts for the differences measured between different observations, different channels and 
different spectral ranges where lines of the species concerned are seen, but observed through 
the same aperture (see the discussion in Hebrard et al. 2002). Note that the total O I or 
N I column density evaluations made with the LWRS+MDRS are in both cases lower than 
those made with the other apert ures. We do not know the cause of this apparent systematic 
effect, but included it in the overall uncertainty on our result. 

The quadratic sum of statistical and systematic uncertainties then leads to: 

log N(0 i) tot = 14.86 (±0.07) 
log iV(N l)tot = 13.87 (±0.07) 


2.2. Measurement of N(D i) 

In the FUSE bandpass, D I is clearly detected at Lyman /?, weakly at Lyman 7, and 
not at Lyman 5. We list in Table 4 the different estimates of the total column density of D I 
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that have been obtained from HST observations of Lyman a, and which we measure using 
the present FUSE data. Here as well, we assumed that only one component is present on 
the line of sight since the 1 > I lines are optically thin. Again we checked the impact of a more 
complex line of sight structure and found it to be negligible, as expected. Vv e carried out 
the same investigations of systematic effects as for O I and N I, using the approach detailed 
above and in Hebrard et cl. (2002). 

There is however one essential difference between D I and O I or N I in terms of sys- 
tematics. Namely the cont inuum to the D I absorption is provided by the blue wing of the 
corresponding H I absorpl ion, while the continuum around the O I and N I is very smooth. 
It is thus necessary to measure the possible systematic tied to the estimate of the stellar 
continuum on the measurement of A'(D i). In order to do so, we have measured A r (D i) with 
and without beforehand correction of the data by a theoretical NLTE stellar continuum (de- 
scribed in Section 3). Some corresponding fits are shown in Fig. 2, and the best-fit to the 
HIRS data is shown in Fig. 4. The effect of any assumption on the stellar continuum was 
found to be negligible as compared to the statistical uncertainty. 

However, we have found that the value of N ( D i) measured if the H I lines are excluded 
from the fit and their bl le wing modeled by a polynomial is higher than when the H I 
are included in the fit. This systematic effect may result from the following: in the latter 
(standard) approach, where both D I and H I are fitted together, the position of the D I 
absorption is tied to that >f H I, i.e., they have the same radial velocity and their zero point 
wavelength offsets are the same as both lines appear in the same sub-spectrum. However, 
in the former approach, i.he velocity of D I is not subject to this constraint. These two 
approaches are denoted in Table 4 as “D I & H i” (standard approach) and “D I, no H I” for 
the other one. One should note that a discrepancy in radial velocity between the D I line 
and its counterpart H I could be attributed to small scale ripples in the FUSE wavelength 
calibration or to the pres- rice of weak H I absorbers which would shift the position of the 
H I line with respect to that of the D I line in which these absorbers would not be felt. In 
this respect the approach in which D I is fitted independently of H I seems more adequate; 
however the interpolation of the blue wing of the H I line by a polynomial is not always 
straightforward. We should mention that this systematic effect is not fully understood and 
does not seem to be always present on different lines of sight. In any case it is included in 
our estimate of the systematic error. 

Finally a small additional systematic effect is found to be attached to the uncertainty 
in the shape of the LSF, i.e. whether it is a single gaussian or a double gaussian with wide 
wings (indicated as “double LSF” in Table 4). As mentioned previously, a single gaussian 
LSF has a free FWHM in the fit, while a double gaussian LSF has free amplitude ratio 
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between both components, free FWHM for each, but separation between both components 
fixed to zero pixel. The impact of the shape of the LSF on the measured N(D i) is shown 
in Figure 5 where the Ay; 2 is plotted as a function of logiV(D i). In practice this curve is 
calculated by fixing log./V(D i) to a given value, finding the best fit x 2 f° r this value with all 
other parameters free, and plotting the difference between this x 2 and the best x 2 obtained 
for all possible values of log/V(D i). In this figure, the curvature of a Ay; 2 curve gives the 
statistical error, as usual, while the relative shift between two curves corresponding to two 
different sets of assumptions is interpreted as an estimate of the lcr systematic error tied to 
the uncertainty on the assumptions. 

We thus conclude that the total D I column density on this line of sight measured using 
the HIRS data is: 


log N(D i)™ RS = 13.39 ± 0.07 ± 0.06 

where the errors are first statistical (2a) and second systematics (95.5% c.l.). 

Combining these errors leads to: 

logiV(D l)™ RS = 13.39 ± 0.09 

As can be seen from the comparison between the evaluations of iV(D i) made for the 
different apertures, there is no significant systematic effect tied to the aperture, contrary to 
the analysis for O I and N I. The average N(D i) measured through the different apertures 
then reads: 


logiV(D l) tot = 13.40 ±0.07 

This value agrees with all previously measured values of N ( D i) shown in Table 4, 
except with that of Sahu et al. (1999), using the STIS-Ech#l data, which gave log./V(D i)= 
13 . 55^0 Us- ^ thus appears that the discrepancy between this value and the others should 
be attributed to the dataset used, in agreement with the conclusion of Vidal- Madjar (2001), 
and the debate around the value of N(D i) is now settled. 


3. The velocity structure of the line of sight and Af(H i) 

In this Section we analyze new high resolution high signal- to- noise STIS data of G191- 
B2B in order to determine the number of components on the line of sight and provide a new 
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estimate of the total neul ral hydrogen column density. As mentioned in Section 1, there 
exists a controversy in the literature about the total number of absorbing components, as 
Sahu et al. (1999) claimed to see two only, whereas Vidal-Madjar et al. (1998) claimed 
that three at least were present. In this Section, we present concrete evidence in favor of 
the latter. Previous measurements of Af(H i) are compiled in Table 5; these values scatter 
around a (non-weighted) mean logiV(H l)= 18.34 ± 0.03 (the dispersion corresponds to 
the non-weighted dispersion of the individual measurements around the mean). Individual 
measurements tend not to agree with each other, and therefore a new 7 independent estimate 
is useful. Moreover w 7 e argue in this Section that the total A(H i) toward G191-B2B is, as 
a matter of fact, much less well known than previously thought, since additional weak hot 
absorbing components may strongly affect the column density estimate from Lyman a. The 
new STIS data analyzed here represent about three times the total exposure time cumulated 
hitherto tow 7 ard G191-B2B. 

All STIS observations were extracted following the method of Howdt and Sembach (2000). 
We selected spectral lines corresponding to N I, O I, Si II, Si ill, C II, S II, S III in the E140H 
echelle configuration and Fe II lines in the E230H one, both at R ~ 90,000 resolving power, 
as reported in Table 6. All these species are seen in one or more spectral lines, strong and/or 
weak, and cover a w 7 ide i ange of atomic masses, so that this data should provide strong 
constraints on the structure of the line of sight as well as on the component to component 
physical state and thermal and non-thermal broadening of each. 

In order to verify that no photospheric line is blended with one of the interstellar lines 
analyzed, w r e have computed a synthetic spectrum of G191-B2B w 7 hich includes all species 
observed in the atmosphere of G191-B2B, i.e., C, N, O, Si, P, S, Fe, and Ni (Lanz et al. 
1996; Vennes et al. 1996). We employ the program TLUSTY developed by Hubeny & Lanz 
(1995) to compute a NLTE metal line-blanketed model by using the atmospheric parameters 
of Barstow et al. (1998), i e., log g = 7.4 and T e g = 54 000K, and the abundances determined 
by Barstow et al. (2001). From the comparison of the observations with the model calculation 
it w 7 as possible to select the lines w 7 here no or very weak photospheric features are coincident 
with the interstellar ones, as for the FUSE data. Furthermore it is possible that some 
photospheric lines are present but not predicted by the model; we thus fitted simultaneously 
as many different spectral lines of each species as available, in order to identify and minimize 
the effect of unpredicted photospheric lines, as w 7 as done for the analysis of the FUSE data. 
This approach also reduces the possible impact of instrumental defects w r hich may be present 
in some area of the detector but not in others. 

Following the general method of analysis described in the previous Section, w T e split the 
STIS data into sub-spectra, each of them centered on one spectral line and having typical 
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width ~ 0.3A. All sub-spectra are then fitted simultaneously. The wavelength zero point 
offsets are left free to vary during the fit. Note that shifts in the wavelength calibration from 
region to region could still be present even after a careful wavelength scale calibration. As 
an example, we detected an erroneous 4km s _1 wavelength shift between the GHRS Si ill 
line position and that in the STIS data. After comparison with GHRS first order grating 
observations, we confirmed that the STIS calibration was correct, and the GHRS calibration 
erroneous. We found that the average relative shift between two spectral regions in the 
STIS data is ~ —0.1 ± 1.1km s -1 (2 a error). An rms error of 0.5km s -1 in the STIS 
calibration is indeed compatible with the ~3.3km s -1 resolution of the instrument. The 
stellar continua in each region are interpolated by low order polynomials (except for Lyman a 
where the theoretical continuum is used in some cases, see below), and the background level 
in each region is determined using the closest saturated line. The typical background flux is 
—2.6 ± 2.4% (2cr error) of the continuum for the E140H data, and —0.10 ± 4.2% (2a error) 
for the E230H data. Finally in some cases (see below for details), the LSF is left free to vary 
and corresponds either to a simple gaussian or to a double gaussian. Again, we emphasize 
that the uncertainty on these instrumental characteristics is contained in the final error bars 
given below, which we determine using a A% 2 method, since the corresponding parameters 
are left free to vary during the fit. 

In the case of the STIS observations made with the E140H and E230H echelle modes, 
we used the tabulated LSF, corresponding to the slit used for these observations (0.2 x 0.2 
arc sec). This is an important issue because these LSFs possess wings that could have 
some impact on the precise determination of the physical parameters we are aiming at. 
To test the quality of the fits we also tried a single gaussian LSF with free FWHM. We 
found excellent fits with this freely variable width single gaussian LSF for both echelle 
datasets. We find for the E140H an average FWHM= 2.91(±0.44) pixel (2 a) and for E230H, 
FWHM= 2.38(±0.47) pixel (2a). These determinations are certainly compatible with the 
tabulated LSF in terms of “average” widths, but are different in terms of shape. They 
correspond to 80,000 and 100,000 resolving power respectively for the E140H and E230H 
spectra. Finally, we also tested double gaussian LSFs but found no significant improvement. 
We thus decided to present results with both types of LSFs (free gaussians and tabulated) 
in order to show the stability and robustness of our column density determinations. As we 
will see, such a change has negligible impact on the derived total H 1 column density. 

Following the analysis of the FUSE data, we evaluate the accuracy of the intrument 
estimate of the individual pixel errors by completing low order polynomial fits over relatively 
smooth and flat continuum sections close to each of the fitted spectral lines. These simple 
fits give a total x 2 =921.37 (sum of all x 2 i n the different portions of continuum analyzed) 
for 661 degrees of freedom, which lies well above the upper 95.5% confidence level limit of 
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723.9 for a y 2 distributioi with 661 d.o.f. Therefore we rescale all x 2 by a common factor 
921.37/600.63 = 1.53, where 600.63 corresponds to the lower 95.5% confidence level limit for 
661 d.o.f. 

Finally, in the course of the analysis, the physical characteristics of each absorbing 
component (velocity v, temperature T and micro-turbulent broadening f) are determined; 
however we will not provide a detailed analysis of these characteristics and their uncertainties 
since we are primarily interested in the total number of absorbing components and the total 
neutral hydrogen column density. 


3.1. Number of absorbing components 

In order to settle t he debate on the total number of absorbing components, we have 
studied two and three component solutions. Simple eye inspection of the data indicates 
that at least two components are present on the line of sight; as mentioned previously, the 
red component can be identified with the LIC, a moderately neutral region, while the blue 
component is clearly more ionized, and, as we argue, more complex than a single absorber 
(see Vidal-Madjar et al, 1998 and Sahu et al. 1999). 

The result of this analysis is summarized in Table 7, ■which gives the y 2 values obtained 
for the fits of various spectral regions, and for all regions (excluding Lyman ct), for two 
and three component solutions, using either the tabulated STIS LSF or a freely varying 
single gaussian LSF. It is clear that the improvement in y 2 when going from two to three 
components is extremely high whichever LSF is used. This is particularly true for Fe II in 
which case the improvement can be seen directly on Figure 6. Indeed if one looks closely at 
the Fe II lines, one can see an asymmetry of the line profile of the bluer component which is 
common to all Fe II lines. This constitutes clear evidence for the presence of three absorbing 
components and not two. One sees this third component in Fe II and not in other lines since 
Fe II is the heaviest species, and its line widths are thus the smallest. Moreover its lines 
have been observed through the echelle presenting the higher resolving power (100,000). 
This suffices to reveal the double structure of component B: components B1 and B2 are 
separated by 3.8 ± 0.5 km s -1 (2cr), i.e. of the order of the E140H spectral resolution but 
more than that of the E230H. The radial velocities of the three components is estimated 
as wbi = 7.7 ± 0.5km s -1 . ub 2 = 11.5 ± 0.5km s _1 and i>lic = 19.4 ± 0.5km s -1 , with la 
errors due to the wavelength calibration. The LIC velocity is in perfect agreement with its 
projected velocity on the line of sight. 

One can also see further evidence for the double structure B1-B2 in the S ill 1190.2.4. 
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line shown in the upper right panel of Fig. 7. This line is seen only in the bluest component 
Bl, which reflects the difference in ionization between the various components; its spectral 
position is quite precisely known as it is located close to the Si II line at 1190.4.4 and is seen 
in two different spectral orders (see Table 6). Note that the position of most lines is actually 
very well controled, which is important to determine the number of absorbing components. 
In effect, the sharp geocoronal lines of N I, 0 I and O I* are in several locations (see Fig. 7), 
notably close to interstellar lines of N I, O I and Si II. Other Si II lines or Fe II lines are 
well constrained since in each of these species the absorption due to the LIC is sharp and 
pronounced, and serves as a useful spectral reference. 

In order to test the hypothesis that the presence of a third component is not required by 
the fit, we perform an F— test. This test uses the Fisher-Snedecor law which describes the 
probability distribution of the x 2 ratio. Here we test the probability that the decrease of the 
X 2 due to the inclusion of a third component could be due to the increase of free parameters 
and not to actual information contained in the data. The result is shown species by species 
and for all species fitted simultaneously in Table 7, and gives a probability < 1CT 4 that this 
third component is not required by the data. The F — tests performed species by species also 
confirm our previous impression that Si II and Fe II are the two most sensitive species to the 
presence of this component. We also note that the need for a third component is present 
whichever LSF is used. 

We note that the x 2 values corresponding to the use of the tabulated STIS LSF are too 
high, as they all lie above the upper 95.5% confidence level limit for statistical fluctuation 
(an except being Fe II for three components). In contrast, the x 2 values for a freely varying 
single gaussian LSF are satisfactory, since only O I and Si n stand out with x 2 beyond the 
95.5% confidence level for the three component solution. Their high x 2 propagates into the 
total x 2 summed over all windows. However the Si II 1190A, 1193A and O I 1302.4. lines 
are saturated, and small errors in the estimates of the error bars at the bottom of the line 
could potentially explain large differences in x 2 - In effect our procedure of rescaling the x 2 
amounts to rescaling all error bars by a common factor, which may not be correct for pixels 
at the bottom of saturated lines for the following reason. Error bars on flux values arise as a 
combination of background noise, fixed pattern noise and Poisson noise. Background noise 
dominates for pixels with flux values close to zero, while for high S/N data Poisson noise or 
fixed pattern noise would dominate in the continuum. If the errors in the estimation of the 
noise array for the background component and for the continuum component (fixed pattern 
or Poisson) are different, then the procedure of rescaling error bars with a common factor 
is, strictly speaking, incorrect. Finally, since this rescaling factor was estimated from the 
continuum, this rescaling should be correct for all pixels whose error is not dominated by the 
other noise component (background). We also note that this 0 I line has a complex profile, 
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with two geocoronal and one nearby photospheric feature (see Figure 7). Nevertheless, we 
decided to keep this line in our analysis as it adds important constraints on the B1 and LIC 
component T and £ parameter evaluations. 

Using the F— test, we can determine the probability that the use of the freely varying 
Gaussian LSF is not required (instead of using the STIS LSF). In all cases, it is found that 
this probability is smaller than 10~ 4 , and therefore we conclude that the gaussian LSF with a 
width that varies from spectral region to spectral region is required to fit the STIS data. This 
result may be unexpected but it can find a simple explanation in the fact that the tabulated 
STIS LSF is an average LSF given over the whole spectral range, while our gaussian LSF 
has a FWHM that varies with wavelength. For standard datasets it is probably sufficient 
to use the STIS LSF but the very high quality of the present data necessitates higher order 
corrections to this LSF. In particular we noted a slight broadening of the LSF near the order 
edges. Furthermore, the wings of the tabulated LSF are not so important and the single 
gaussian LSFs are sufficient to properly fit the data. 

We thus conclude the t three absorbers are present on the line of sight to G191-B2B. 


3.2. Lyman a 

The background level is directly measured at the bottom of the strongly saturated 
Lyman a line. It is found to be of the order of — 1.0(±1.0)% (2a) of the nearby continuum, i.e. 
within a few percent. Variations in the photospheric continuum over the interstellar Lyman a 
line presents an additional difficulty. Indeed the instrument sensitivity over this spectral 
region is evaluated by using model spectra of white dwarfs and comparing them to the data, 
and one of the calibration WDs used is G191-B2B itself. Therefore it is difficult to separate 
the instrumental effects from the photospheric intrinsic profile in the spectrum processed 
through the STIS pipeline. For this reason, we follow the approach detailed in Vidal-Madjar 
et al. (1998) and further discussed by Vidal-Madjar (2000), and evaluate simultaneously all 
parameters and the stellar continuum, which we model in two different ways. In the first 
approach, we model this continuum as a polynomial whose parameters are adjusted during 
the fitting procedure, and which represents the real photospheric profile times the instrument 
sensitivity; we denote this fitting model “ IF , which stands for uncorrected (as the data have 
not been corrected prior to the fit by a stellar model). In the second approach, we correct 
the data prior to the fit be a stellar model and keep a free low-order polynomial during the 
fit which models the inaccuracy of the stellar model and the variation of the instrument 
sensitivity; we denote this approach by “C”. The calculated stellar model was shifted by 
24.56 km s _1 in order to be consistent with the velocity of the photospheric features clearly 
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detected in N v and Si in (see Figure 7). These two approaches have been followed in parallel 
to estimate the uncertainty in the modeling of the stellar continuum on the final solution. 

We use two models to compute the photospheric Lyman a line profile. First, we use 
the program SYNSPEC (I. Hubenv 2000, private communication) in conjunction with the 
TLUSTY NLTE metal line-blanketed model described above. This version of SYNSPEC 
contains the Lemke’s Stark broadening tables for hydrogen, which were computed by Lemke 
(1997) within the framework of Vidal, Cooper, & Smith (1973). Second, we use Detlev 
Koester’s LTE code and his best fit parameters derived under the assumption of pure-H 
LTE (Te ff = 60880 K and log g = 7.59). Figure 8 shows the NLTE and LTE Lyman a line 
profiles and illustrates that NLTE effects are significant only in the core of the line, which 
is formed high in the atmosphere where departures from LTE are important (Wesemael et 
al. 1980). However, as is shown in Fig. 8 the observations do not contain any information 
in that central region of the photospheric line simply because it is lost at the bottom of the 
saturated interstellar Lyman a line. Our fitting procedure will however be able to test the 
slightly different slopes in the wings of the models. 

An additional correction made by a low-order polynomial is needed not only to take into 
account uncertainties in the instrument sensitivity but also because when one fits a section 
of the G191-B2B stellar continuum devoid of stellar or interstellar absorption, it is necessary 
to model the continuum with a polynomial of degree of order 1 to 3 for a spectral width 
of order 0.3 A (see for instance Fig. 7). For a spectral region of width ~ 3A as is the case 
for Lyman a , a 6 th order polynomial provides a satisfactory approximation to the stellar 
continuum. 

The best-fit solutions obtained for 3 components on the line of sight for all lines in 
the STIS domain including Lyman a. and for various models (prior correction or not of 
the Lyman a stellar continuum, type of LSF) are summarized in Table 8. The zero point 
wavelength of the Lyman a spectral domain was found to be compatible with the other 
regions well within the 1 a uncertainty of 0.5km s -1 . 

Concerning the fit with a free double gaussian (model 5 in Table 8), we found that the 
wings of the broader gaussian (about 5 to 10 times larger than the narrow one) contribute in 
amplitude (relatively to the amplitude of the narrow core gaussian) to less than 1% in all but 
5 spectral windows, less than 2% in 3 of these 5 and to about 3% in the 2 remaining ones. 
Since this effect is marginal, we kept only as comparison tests the use of both the tabulated 
STIS LSF or the simple gaussian LSF. All these effects u r ere included in the evaluation of 
the statistical errors. 

It appears that different values of iV(H i) are obtained depending on the model used; in 
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particular the comparison of models 1 and 2 on the one hand, and of models 3 and 4 on the 
other hand shows that the value of A T (H i) derived depends whether the stellar continuum at 
Lyman a has been correcled or not beforehand by a theoretical stellar profile. One obtains 
logJV(H i)=18.37 for an ur corrected profile, and logiV(H i)=18.32 for a corrected stellar pro- 
file. Since there is no significant difference in x 2 between these solutions, this discrepancy 
must be interpreted as a systematic uncertainty tied to the choice of modeling of the con- 
tinuum. A Ay 2 analysis around each of these solutions give a statistical error of ±0.01 dex 
(lcr), so that the value of ,V(H i) is dominated by the above systematic error. Furthermore 
after correction of the stellar continuum by an LTE calculation, which as shown in Fig. 8 
represents an intermediab' situation between NLTE and non-corrected stellar profiles, the 
value of N ( H i) derived was found to be an intermediate value between the above two. 

However a further systematic error appears if additional weak hot components are 
present on the line of sight., and its investigation is the subject of the following Section. 


3.3. Systematic uncertainties and additional hot components 

One cannot exclude 1 he presence of weak H I components which could perturb or bias 
the measurement of 1V(H i), but could not be detected in any other species due to their 
weak column density. Su ch absorbers could arise as the signature of high velocity shocks 
(Cowie et al. 1979), or as cloud interfaces with the hot gas within the local ISM (Bertin et 
al. 1995) or as “hydrogen walls”, i.e., the shock interaction between the solar wind (or stellar 
wind) and the surrounding ISM (Linsky 1998). The latter possibility has been modeled by 
Wood et al. (2000), and a prediction for the line of sight toward G191-B2B is shown in 
Fig. 9 overlaid on the STIS data. It shows that most of the expected absorption should take 
place in the saturated core of the interstellar line, but some weak absorption (~ 5%) may be 
present, extending over several tenths of an Angstrom on the red side of the line, due to the 
neutral hydrogen atoms seen behind the shock in the downwind direction w r here G191-B2B 
is located. 

We have thus investigated the possible impact of additional weak hot components on 
the determination of N(H i). In order to do so, we have added one or tw T o additional 
components in H I only, and performed the fit of all lines as before with an NLTE correction 
of the Lyman a profile and using a free single Gaussian LSF (model 3). In order to constrain 
the presence of these exlra absorbers, w r e have added to the fit the four Lyman f3 lines 
observed through the FUSE HIRS aperture. For each new model characterized by the 
number of additional components, we have measured iV(H i) using a Ay 2 analysis. The best 
fit solutions with zero, one and t.wo extra absorbers are shown for Lyman a and Lyman 8 in 
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Fig. 10. We have found that this introduction of additional components leads to a substantial 
decrease of the total N( H i) together with a significant improvement of x 2 - In particular, 
we found that the best fit solution for one extra component leads to Ax 2 = 26 (1961 d.o.f.) 
and logiV(H l)=18.24, and two extra components lead to Ax 2 = 39.4 (1958 d.o.f.) and 
logA r (H i)=18.11. The x 2 values quoted have been rescaled by a factor 1.53 corresponding 
to the STIS data obtained in the previous section, and which remains close (20%) to the 
factor measured in the FUSE data range. Moreover, if a Ax 2 analysis is performed for each 
model, the curve obtained flattens for low values of N( H i) down to logiV(H l)=18.0; in other 
words, the space of solutions becomes degenerate in the low N( H i) region, and consequently 
the value of N( H i) measured from the profile fitting of the Lyman a region is subject to 
a large systematic uncertainty. This behaviour can be explained as follows. In the three 
component best fit solution, most of the hydrogen is contained in the LIC whose absorption 
makes the red wing of Lyman a. In the solutions with extra absorbers, the combined profile 
of these hot components can replace the damping wing contribution of the LIC, as shown in 
Fig. 10. 

In order to scan the parameter space, the iV(H i)lic h as been forced to take various 
values, since the profile fitting code cannot force the total H I column density to take a 
particular value. When the LIC iV(H i) becomes negligible with regards to the total iV(H i), 
we force the next dominant component to take various values of JV( H i). In this way it is 
possible to scan the parameter space of the total iV(H i). One finds that once the LIC has 
become negligible, component B2 makes most of the total N(H i), and that when the total 
iV(H i) is taken below 18.0, the x 2 rapidly increases, which gives the final error bar on the 
total JV(H i) estimate. 

At this stage, one should point out that these best-fit solutions are quite difficult to 
find. The impact of additional weak absorbers on the profile fitting of N(R i) had already 
been studied by Vidal-Madjar et al. (1998) but no significant effect had been found, as the 
above solutions had eluded detection. The above solutions have actually been obtained only 
in a late stage of the present work. We interpret this as evidence for the fact that the x 2 
surface becomes complex when additional unconstrained absorbers are introduced in the fit 
of the fully blended Lyman a profile and optimization is then delicate. 

One should also remark that we have modeled these extra absorbers using Voigt profiles, 
which implicitly assumes a Maxwellian distribution for the velocities of the atoms. However 
if these absorbers correspond to interstellar structures similar to hydrogen walls, the latter 
assumption is incorrect and the overall profile should be closer to that modeled by Wood et 
al. (2000) and shown in Fig. 9. For this reason, the above modeling and solutions should be 
interpreted as a first approximation. For this reason as well, we do not discard the values 
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for iV(H i) obtained in th- previous section without additional absorbers, even though the 
presence of extra absorbers induce a significant gain in x 2 - At the worst the final value for 
/V(H i) will be more conse rvative. With these remarks in mind, we conclude that the total 
hydrogen content toward < 1191-B2B should be contained in the interval: 


18.00 < log A r (H i) < 18.37. 

In the absence of any indication on the distribution of the errors for A r (H i), we consider 
the above two extreme values as 2cr limits since they correspond to rather extreme solutions, 
and quote the following value for N(H i): log./V(H l)= 18.18 ±0.09 (ler error). One should 
note that this effect of additional absorbers had not been found in previous studies of the 
total hydrogen content toward G191-B2B. Therefore it is important to remark that the above 
large uncertainty must also affect the previous values of AT(H i). In other words, the above 
results supersedes any previous estimate of N(¥L i) made toward G191-B2B using profile 
fitting of the Lyman a line. One may also wonder whether extra absorbers may affect the 
determination of AT(H i) toward other stars; this important question is the subject of a 
forthcoming paper (Vidal Madjar &; Ferlet 2002). On the other hand, the above value for 
N(H i) should be contrasted with the value derived from EUVE observations and modeling 
of the atmosphere of G191-B2B which gave logAf(H i)=18. 315 ± 0.013 (Dupuis et al. 1995, 
2 a error), logA r (H l) = 18.32 (Lanz et al. 1996, no error bar quoted) and the recent detailed 
measurement logAf(H l) = 18.30 ± 0.09 (Barstow & Hubeny 1998; this value also agrees with 
the more recent work of Barstow, Hubeny fc Holberg 1999). These values would tend to 
indicate that the contamination of the Lyman a profile by weak hot absorbers (if any) is not 
important, as they agree w ith the values we obtained without including such extra absorbers. 
However, as noted in Barstow & Hubeny (1998), these fits of the EUVE spectrum of G191- 
B2B typically lead to very large and unexplained reduced x 2 > implying that the overall fit 
is not yet satisfactory in spite of the use of sophisticated W D atmosphere models, and that 
some unknown effects have yet to be accounted for. 


4. Discussion and conclusions 

We have measured the total column densities of D I, N I and O I toward G191-B2B 
using unsaturated absorption lines of these elements in high quality FUSE spectra. After a 
careful examination of the possible systematic uncertainties tied to the choice of the stellar 
continuum and to the instrumental configuration, we have derived the following column 
densities with 2a uncertainties: 
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log N(D i) tot = 13.40 ±0.07, 
log N(0 i) tot = 14.86 ±0.07, 
log iV(N i) tot = 13.87 ±0.07. 


We have also analyzed new high signal-to-noise ratio high resolution STIS observations 
of G191-B2B and provided concrete evidence for the presence of at least three interstellar 
absorbing components on the line of sight by analyzing the interstellar absorption lines of 
N I, O I, Si II, Si in, S ill and Fe II present in the STIS bandpass. We have also measured 
the total hydrogen column density on the line of sight using the velocity structure derived 
from the above metals. We have performed an exhaustive study of systematic effects on 
the value of N(E i). In particular we have discovered a new major source of uncertainty 
on iV(H i) tied to the possible presence of additional weak hot absorbers whose combined 
absorption profile can contribute significantly to the wings of the blended Lyman a profile. 
The column density of these absorbers is small compared to the other main components, 
and they would not be detected in any other species than H I, but their contribution to the 
Lyman a absorption profile can reduce significantly the total H I column density measured 
from the profile fitting. In order to constrain their impact, we have analyzed simultaneously 
Lyman a and the higher order Lyman lines, and concluded that the best value of N ( H i) 
toward G191-B2B is: 


logiV(H l) tot = 18.18 ± 0.18 (2(7 error) 

We emphasize that this uncertainty is a systematic uncertainty which had gone unno- 
ticed before. Therefore the above result supersedes previous estimates of JV(H i) toward 
G191-B2B obtained from the profile fitting of Lyman a. A detailed analysis of this un- 
certainty and its consequences on iV(H i) determinations toward other stars is discussed 
in a companion paper (Vidal-Madjar k Ferlet 2002). We thus derive the following neutral 
abundance ratios toward G191-B2B, with 2a uncertainties: 


(D/H) tot 

= 1.66±g;| 10~ 5 

(D/0) tot 

= 3.49 ± 0.78 10~ 2 

(D/N) tot 

= 3.41 ± 0.76 10 _1 

(0/H) tot 

= 4.79tfj 10~ 4 

(N/H) tot 

= 4. 90^i g 10 -5 


( 1 ) 
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Most of the uncertainty in the above result results from the systematic uncertainty on 
the JV(H i) determination This clearly shows the importance of measuring accurate (D/O) 
and (D/N) ratios in the interstellar medium instead of abundances relative to hydrogen, as 
emphasized by Timmes e' al. (1997). Interestingly if one uses the recent measurement of 
jV(H i) from the modeling of the atmosphere of G191-B2B and the fit of the EUV E spectrum, 
loglV(H i)= 18.30 ±0.09 (2(j), one finds (D/H) tot = 1.26l^ x 10" 5 (2a). At this stage, 
however, due to the uncertainty inherent to the modeling of the white dwarf atmosphere, it 
is probably more conservative to use the interstellar determination for A^(H i), and therefore 
the previous value of the D/H) ratio. 

The above new value for the (D/H) ratio agree with the range of values measured by 
Linsky (1998) toward a dozen stars of the LISM and with the values previously derived 
toward G191-B2B. However the discrepancy between previous estimates of the (D/H) ratio 
toward G191-B2B and the LISM average D/H ratio has disappeared due to a revision of 
the uncertainty on the esl inflation of the total H I content. A detailed interpretation of this 
(D/H) value and of the accompanying (D/O) and (D/N) ratios and their implications is 
provided in a companion oaper by Moos et al. (2002). 
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971.5 972 972.5 973 


Wavelength (A) 

Fig. 2. — FUSE MDRS spectrum of G191-B2B around Lyman 7 (segment SiC2A). The 
solid line in each panel shows the best fit obtained by a simultaneous profile fitting of 
all H 1, D I, O I and N I lines of all FUSE datasets. Clearly these observations do not 
resolve the multiple absorbers along the line of sight. In the upper panel the continuum is 
interpolated by a smooth polynomial, and three interstellar components corresponding to the 
velocity structure of the line of sight derived from the higher resolution STIS observations (see 
Section 3) were considered; in the middle panel, the stellar continuum is also interpolated 
by a polynomial, while only one absorbing component is considered. In the lower panel 
the stellar continuum has been corrected by a theoretical NLTE stellar profile, the residual 
continuum is modeled by a polynomial and only one absorbing component has been assumed. 
See text for details. 
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Wavelength (A) 

Fig. 3 — FUSE spectra of G191-B2B around the N I multiplet at 954A. The aperture, 
segment and resolving power R corresponding to each dataset are indicated on the figure. 
The line at 954.1 A was excluded from the profile fitting since it is likely contaminated by 
a photospheric feature, as indicated by theoretical modeling of the stellar continuum, and 
as shown clearly in the upper and middle panels; in the lower panel, the apparent absence 
of this line may result from a statistical fluctuation or from a detector feature. The higher 
resolution of the HIRS aperture data is apparent. 
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Radial velocity (km/s) 


Fig. 4 . — FUSE HIRS spectra of Lyman f3 and Lyman 7, recorded on the segments indicated 
in the figure. The solid line shows the final best-fit solution to all lines of H I, D I, N I and 
O I in all HIRS data fitted simultaneously. The weak emission line at the bottom of the 
Lyman lines is due to H I geocoronal emission; the pixels affected by this emission are not 
considered in the fit. The stellar continuum is interpolated by a smooth polynomial. See 
text for details. 
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Fig. 5. — Curve of Ay' 2 deviations around the best fit \ 2 as a function of logiV(D i). Only 
the curve for the HIRS data is shown here; statistical error bars for other aperture data 
were obtained using the same method. The A% 2 values shown here have been rescaled to 
compensate for uncertainties in the individual pixel errors (see text). The various curves 
correspond to different models for the fit: in thick line, the D I and H I are fitted simultane- 
ously with a LSF modeled as a single gaussian with free FWHM; in intermediate thickness, 
the same approach but with a double gaussian LSF, with free amplitude ratio and FWHMs; 
in thin line, H I is excluded from the fit, the continuum to the D I absorption is modeled by 
a polynomial, and the LSF is a simple gaussian. The curvature of these curves give the the 
statistical errors while their relative shifts give an estimate of the overall systematic uncer- 
tainty related to the different models. The hatched area shows the final 95.5% confidence 
level error for logiV(D i) measured using the HIRS data. 
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Fig. 6. — The STIS data covering the three Fe II lines observed with the E230H Echelle 
grating. All fits shown here were made using a freely varying single gaussian LSF ; the Fe II 
lines were fitted simultaneously with all other species including Lyman a (see text). The left 
panels show the fits with two absorbing components only, and the right panels show the fits 
using three components. One can clearly see the same asymmetry of the bluer component 
B in all Fe II lines which reveals the presence of complex substructure in this component, 
hence the need for 3 absorbers in total. The x 2 corresponding to these fits is 227.9/139 for 
the two component solution, and 123.3/136 for the three component solution; the gain is 
clear and comforts the visual impression. The error bars on each pixel are of order of 2% of 
the continuum level. 
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Fig. 7. — The STIS data with the best-fit three component solution using the tabulated 
STIS LSF; only 8 sub-spectra out of 19 in total are shown here, and are labeled with the 
central spectral line wave length and Echelle order: the N I 1200A triplet along with the 
nearby geocoronal absorptions (marked “Geo.”, in two different locations because two data 
sets taken at two different epochs are here averaged), the O I line with the corresponding 
geocoronal O I absorptions and a nearby photospheric feature (noted “Phot.”), three of the 
Si II lines with either the nearby S III feature or the O T geocoronal absorption and the Si III 
line with an additional broad Si III photospheric line located at the same velocity shift as 
the N v photospheric line (which was fitted simultaneously). Note also that the S III lines 
clearly detected in two different orders near 1190.2A are also spectrally well located. The 
fact that they show up only in the bluest component B1 is an additional reason why 3 and 
not 2 components are needed along that line of sight. 
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Wavelength (A) 


Fig. 8. — The STIS data over the Lyman a spectral region along with the LTE (dashed 
line) and NLTE (dotted lines) photospheric profiles. Three different NLTE profiles for at- 
mospheric parameters in the ranges 54000 K < T e ff < 55000 K and 7.5 < log <7 <7.6 are 
shown. The LTE calculation uses T e ff = 60880 K and log g = 7.59. The radial velocity for 
the photospheric profiles is set to +24.56 km s -1 . The two thin vertical lines delimit the 
wavelength range in which no information on the stellar profile is contained in the data. The 
fitting procedure thus tests only the difference between the wings of the models and not their 
cores whose difference is more pronounced. 
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Fig. 9. — The Lyman a interstellar absorbtion observed with STIS with the hydrogen wall 
model of Wood et al. (2000) for the line of sight of G191-B2B overplotted in dashed line. 
The STIS data has been normalized by the best-fit continuum. 
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Wavelength (A) 

Fig. 10. — Best-fit solutions for all lines in the STIS data plus four Lyman fi from the FUSE 
HIRS dataset using three main interstellar components plus zero, one and two extra H I 
components (upper, middle, and lower panels respectively). The total absorption profile 
is shown as the solid line, and the dotted and dashed lines show the individual profiles 
(convolved with the LSF); the thick dashed lines in particular show the contribution of the 
extra absorbers. The solutions are labeled with the total rescaled x 2 (for 1964, 1961, and 
1958 d.o.f. respectively for the upper, middle and lower panels) and total H I column density 
(in log). 
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Table 1. Strongest N I and O I spectral lines in the FUSE domain. 


Wavelength i A) 

Element 

f 

A (s J ) 

Comment 

1025.4 440 

D i 

0.264 10' 1 

0.190 10 9 


1025.4429 

D i 

0.527 10' 1 

0.190 10 9 


972.2725 

D I 

0.967 10~ 2 

0.813 10 8 


972.2721 

D i 

0.193 10- 1 

0.813 10 8 


949.4848 

D i 

0.465 10~ 2 

0.421 10 8 

Weak 

949.4846 

D I 

0.929 10” 2 

0.421 10 8 

Weak 

1134.9803 

N i 

0.435 10- 1 

0.150 10 9 


1134.4149 

N i 

0.297 lO" 1 

0.154 10 9 


1134.1653 

N i 

0.152 10- 1 

0.158 10 9 


954.1042 

N i 

0.676 10~ 2 

0.330 10 8 

Blend 

953.9699 

N i 

0.348 10- 1 

0.170 10 9 


953. ('549 

N i 

0.250 10- 1 

0.183 10 9 


953.4152 

N i 

0.132 10" 1 

0.193 10 9 


1039.2303 

0 i 

0.920 10" 2 

0.947 10 8 


1026.4757 

0 i 

0.246 10~ 2 

0.111 10 8 

Weak 

1026.4744 

0 i 

0.187 10“ 3 

0.118 10 7 

Weak 

1025.7633 

0 i 

0.201 10“ 3 

0.212 10 7 

Blend 

1025.7626 

0 i 

0.302 10~ 2 

0.191 10 8 

Blend 

1025.7616 

0 i 

o.i69 nr 1 

0.765 10 8 

Blend 

988.7734 

0 i 

0.465 10- 1 

0.226 10 9 

Strong 

988. (-549 

0 i 

0.830 10~ 2 

0.566 10 8 

Blend 

988.5778 

0 T 

0.553 10- 3 

0.629 10 7 

W r eak 

976.4481 

0 i 

0.331 10- 2 

0.386 10 8 


971.7382 

0 i 

0.116 10- 1 

0.585 10 8 

Strong 

971.7376 

0 i 

0.207 10- 2 

0.146 10 8 

Blend 

950.8846 

0 i 

0.158 10~ 2 

0.194 10 8 


948. (>855 

0 i 

0.631 10- 2 

0.100 10 9 


936.6295 

0 i 

0.365 10~ 2 

0.100 10 9 


929.5168 

0 i 

0.229 10~ 2 

0.100 10 9 


925.4460 

0 i 

0.354 10~ 3 

0.459 10 7 

Weak 

924.9500 

0 i 

0.154 10~ 2 

0.100 10 9 
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Table 1 — Continued 


Wavelength (A) 

Element 

/ A (s A Comment 

921.8570 

0 i 

0.100 10- 2 0.562 10 7 Weak 


Table 2. FUSE observation log of G191-B2B. 


Dataset 

Aperture 

-^exp (ksec) 

^exp 

Date 

S3070101 

LWRS 

15.5 

32 

2000.01.14 

P1041202 

MDRS 

15.5 

21 

2000.01.13 

P1041201 

HIRS 

15.5 

32 

2000.11.06 


Table 3. Published and FUSE O I and N I column densities. 


Spectrograph / 
Aperture 

log N(0 i) to t 
2cr error 

logiV(N l) tot 
2a error 

Reference 

HST GHRS-Ech. 

14.84±0.04 

13.90i0.02 

Vidal-Madjar et al. 1998 a 

FUSE LWRS+MDRS 

14.79i0.04 

13.82i0.07 

This work b 

FUSE MDRS 

14.84i0.08 

13.89i0.06 

This work b 

FUSE HIRS 

14.88±0.06 

13.84i0.07 

This work b 

FUSE-AW 

14.86±0.07 

13.87i0.07 

This work b 


a Three components were considered in the fit 
b One interstellar component was considered in the fit 
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Table 4 Published and FUSE D I column densities. 


Spectrograph / 

log A r (D i ) tot 

Reference 

Spect./Aper. 

2a error 

/Model 

HST GHRS-Ech. 

13.43T0.02 

Vidal-Madjar et al. 1998 a 

HST STIS-Ech.#l 

1 3 ^+ 0 - 07 
10 ‘ JtJ -0.08 

Sahu et al. 1999 b 

HST GHRS-Ech. 

13.40i0.04 

Sahu et al. 1999 b 

FUSE LWRS+MDRS 

1 3 41 + 0 - 12 

D I, no H I, this work c 

FUSE MDRS 

13.41i0.09 

D I & H I, this work d 

FUSE HIRS 

13.36i0.08 

Dl&Hi (Fig. 4), this work d 

FUSE HIRS 

13.38i0.06 

D I & H I, double LSF, this work e 

FUSE HIRS 

13.42i0.08 

D I, no H I, this work c 

FUSE All 

13.40i0.07 

this work f 


a 3 free absorbing components assumed in the profile fitting, and stellar 
continuum modeled with a low order polynomial free to vary during the fit 

b 2 free absorbing components in the profile fitting, and stellar continuum 
fixed and modeled by NLTE calculations 

C 1 free absorbing component, stellar continuum modeled by a freely varying 
low order polynomial, and D I fitted alone without the H I line 

d l free absorbing component, stellar continuum modeled by a freely varying 
low order polynomial, and D I fitted simultaneously with the H I line 

e l free absorbing component, stellar continuum modeled by a freely varying 
low order polynomial, and D I fitted simultaneously with the H I line, LSF 
profile modeled with a double Gaussian 


Combination of FUSE models c, d and e 
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Table 5. Published H I column densities. 


log N ( H i) tot (2a) 

Spectrograph 

# Comp. 

Continuum 

Reference 

18.315T0.013 

EUVE 

— 

— 

Dupuis et al. 1995 

18.32 - 

EUVE 

— 

— 

Lanz et al. 1996 

18.36i0.04 

GHRS-G160M 

3 

free 

Lemoine et al 1996 

18.38i0.02 

GHRS-Ech. 

3 

free 

Vidal-Madjar et al. 1998 

18.30±0.09 

EUVE 

— 

— 

Barstow & Hubeny 1998 

18.31i0.03 

STIS-Ech. #1 

2 

fixed 

Sahu et al. 1999 

18.34i0.02 

GHRS-Ech. 

2 

fixed 

Sahu et al. 1999 


Table 6. Spectral lines in the STIS domain used for the line of sight structure study. 

Wavelength (A) 

Species 

/ 

A (s“>) 

Ech. Order 

1334.5320 

C ii 

0.128 

0.288 10 9 

316 

1200.7098 

N I 

0.0430 

0.398 10 9 

350, 351 

1200.2233 

N i 

0.0862 

0.399 10 9 

350, 351 

1199.5496 

N i 

0.130 

0.401 10 9 

351 

1242.8040 

N v 

0.500 

0.336 10 9 

339 

1302.1685 

O I 

0.0519 

0.340 10 9 

323 

1304.3702 

Si ii 

0.0917 

0.107 10 10 

323 

1193.2897 

Si ii 

0.585 

0.409 10 10 

352, 353 

1190.4158 

Si II 

0.293 

0.410 10 10 

353, 354 

1206.5000 

Si III 

1.67 

0.255 10 10 

349 

1259.5190 

S ii 

0.0166 

0.465 10 8 

334 

1190.2030 

S m 

0.0231 

0.651 10 8 

353, 354 

2382.7651 

Fe ii 

0.320 

0.313 10 9 

324 

2586.6499 

Fe ii 

0.0691 

0.272 10 9 

298, 299 

2600.1729 

Fe ii 

0.239 

0.270 10 9 

297 
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Table 7. x 2 comparison for the two and three component solutions using either the 



tabulated STIS LSF or 

freely varying single gaussian LSF. 


Species 

X 2 / d.o.f. 

2 comp. 
STIS LSF 

a x 2 / d.o.f. a 
3 comp. 
STIS LSF 

F— test prob. 
2 vs 3 comp. 

b x 2 / d.o.f. 

2 comp, 
free Gaussian 

a x 2 / d.o.f. a 
3 comp, 
free Gaussian 

F— test prob. 
2 ns 3 comp. 

N I 

395.0/314 

382.3/311 

1.7% 

323.2/309 

309.1/306 

0.3% 

0 i 

214.3/61 

193.1/58 

11% 

113.9/60 

101.4/57 

8.3% 

Si 11 

772.0/468 

651.4/465 

< 0.01% 

678.6/463 

553.6/460 

< 0.01% 

Si m c 

210.6/152 

L92. 8/149 

0.4% 

185.2/150 

170.8/147 

0.75% 

Fe II 

277.6/142 

141.8/139 

< 0.01% 

227.9/139 

123.3/136 

< 0.01% 

All d 

2677.8/1270 

1908.8/1259 

< 0.01% 

2076.7/1250 

1442.5/1239 

< 0.01% 


a the x 2 has been rescal'd to compensate for the inaccuracy of the noise array (see text) 
b The F — test gives the probability that a third absorbing component is not required by the model 
c Fitted simultaneously with the N v line to control the photospheric Si hi line 
d Includes C n, S n and S hi 



40 - 


Table 8. Three component best-fit solutions for all lines in the STIS region including 

Lyman a. 


Fit 

model 

VBl 

(km s _1 ) 

^B2 

(km s -1 ) 

I’LIC 

(km s -1 ) 

logiV (H i) tot 

x 2 

1753 d.o.f. 

^a,c 

7.7 

11.2 

19.4 

18.32 

2469.9 

2 b,c 

7.7 

11.3 

19.4 

18.37 

2470.1 

3 a,d 

7.7 

11.7 

19.4 

18.32 

2006.8 

4M 

7.7 

11.6 

19.4 

18.37 

2001.8 

5 a,e 

7.7 

11.6 

19.4 

18.33 

1987.3 


a data normalized beforehand by NLTE stellar continuum, and stellar 
continuum residuals modeled during the fit by a 6 th order polynomial 

b unnormalized data, with stellar continuum modeled during the fit by 
a 6 th order polynomial 

C LSF corresponds to the tabulated STIS LSF 

d LSF profile modeled with a single Gaussian (free to vary during the fit) 

e LSF profile modeled with a double Gaussian (free to vary during the 
fit) 


